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Standard big-bang nucleosynthesis (BBN) predicts only a trace abundance of lithium and no 
heavier elements, but some alternatives predict a nonzero primordial metaUicity. Here we explore 
whether CMB measurements may set useful constraints to the primordial metaUicity and/or whether 
the standard CMB calculations are robust, within the tolerance of forthcoming CMB maps, to the 
possibility of primordial metals. Metals would affect the recombination history (and thus CMB 
power spectra) in three ways: (1) Lya photons can be removed (and recombination thus acceler- 
ated) by photoionizing metals. (2) The Bowen resonance-fluorescence mechanism may degrade Ly/3 
photons and thus enhance the Ly/3 escape probability and speed up recombination. (3) Metals could 
affect the low-redshift tail of the CMB visibility function by providing additional free electrons. The 
last two of these provide the strongest CMB signal. However, the effects are detectable in the Planck 
satellite only if the primordial metal abundance is at least a few hundredths of solar for (2) and a 
few tenths of solar for (3). We thus conclude that Planck will not be able to improve upon current 
constraints to primordial metaUicity, at the level of a thousandth of solar, from the Lyman-a forest 
and ultra-metal-poor halo stars, and that the CMB power-spectrum predictions for Planck suffer 
no uncertainty arising from the possibility that there may be primordial metals. 

PACS numbers: 98.80.Ft, 26.35. -|-c, 95.30.Jx 



I. INTRODUCTION 

Big-bang nucleosynthesis (BBN) is one of the pillars 
of the hot standard cosmological model. Comparison 
of BBN theoretical predictions to observed abundances 
of the lightest nuclei (D, "^He, '^He and ^Li) uniquely 
determines the only free parameter of standard BBN, 
the baryon-to-photon ratio rj = (5.7 ± 0.3) x lO^^'^, 
or equivalently, given the cosmic microwave background 
(CMB) temperature today Tq — 2.73 K, the baryon 
abundance ^l^/i^ = 0.021 ± 0.001 (see, e.g., Ref. [1]). 
The latest results from CMB anisotropy measurements 
by the WMAP satellite are in excellent agreement, with 
^bh^ = 0.02249t°;°°°^? 0. 

In the standard BBN scenario, elements heavier than 
lithium are only produced with trace abundances F3l. It 
has been shown, however, that significant amounts of 
heavy elements may be produced in inhomogeneous BBN 
(IBBN) models (1H5] . IBBN may take place if some non- 
standard mechanism leads to large baryon-abundance in- 
homogeneities on small scales, which are allowed by cur- 
rent observations. It is possible to adjust the IBBN pa- 
rameters to reproduce the observed abundances of light 
elements, while producing heavier elements with abun- 
dances as large as those in the Sun [51 [7] . More gener- 
ally, it cannot be excluded that some unknown processes 
may lead to a significant production of metals heavier 
than lithium. It may therefore be profitable to inquire 
what empirical constraints to primordial metals might be 
possible. 

Standard methods to constrain metal abundances at 
high redshifts rely on line emission or absorption mea- 
surements, and therefore require some sources to have 
already formed (typically, high-redshift quasars) . The in- 
tergalactic medium (IGM) is therefore already partially 



enriched by metals produced in the first stars, and ex- 
tracting limits on the primordial abundances relies on 
understanding the complex physics of galactic outflows 
and gas mixing and correctly modeling the ambient ra- 
diation field. It would be of great interest to be able to 
probe the abundance of metals before the formation of 
Population III stars that enriched the IGM. A few ideas 
were put forward to probe the metaUicity during the dark 
ages: Ref. [8j suggested using resonant scattering of CMB 
photons off neutral lithium atoms (later shown to be un- 
observable because lithium is kept ionized by redshifted 
Lyman-a photons emitted during primordial hydrogen 
recombination [5]); Ref. [TU] studied the effect of fine- 
structure transitions of heavy elements in atomic or ion- 
ized states on CMB anisotropics; Ref. [TT] considered the 
spectral signatures of carbon and oxygen. In this paper 
we assess whether heavy elements present during primor- 
dial recombination could be detectable from upcoming 
CMB experiments. 

Primordial recombination has recently been the sub- 
ject of a renewed interest, due to the impact of uncertain- 
ties in the standard theory on the predicted CMB tem- 
perature and polarization anisotropy power spectrum. 
Errors in the free-electron fraction Xe{z) as small as a few 
tenths of a percent near the peak of the visibility function 
at redshifts z ~ 1100 would induce biases of several stan- 
dard deviations for cosmological parameters estimated 
from Planck data [121 113] . This accuracy requirement 
has motivated abundant work on radiative transfer in 
the Lyman lines, in particular Lyman-a (see for example 
Refs. [T4HT7] and references therein). The tails of the 
visibility function are less important, but an accuracy 
of about a percent is still needed, which required imple- 
menting an accurate multi-level-atom formulation of the 
recombination problem [T8H2T]. Such a high sensitivity 
to the recombination history can be turned into an asset 
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and serve to probe unusual physics taking place during 
the recombination history as, for example, the presence 
of primordial heavy elements. In this paper, we explore 
this idea, and quantify the impact of neutral metals on 
the Lyck and Ly/3 net decay rates, and of ionized metals 
on the low-redshift tail of the visibility function. 

Below we consider three effects of metals on the recom- 
bination history and thus on the CMB visibility function: 
(1) The removal of Lyman-a photons (and thus accel- 
eration of recombination) by photoionization of metals 
(Section |ll]); (2) the degradation of Lyman-/? photons 
(and thus acceleration of recombination) by the Bowen 
resonance-fluorescence mechanism (Section III); and (3) 



the contribution to the free-electron abundance at late 



times by low-ionization metals (Section IV I. We find 
that effects (2) and (3) provide the biggest impact on 
CMB power spectra. However, the effects are visible in 
Planck only if the primordial metal abundance is at least 
a few hundredths of solar for (2) and a few tenths of so- 
lar for (3). Given that the Lyman-alpha forest [22] and 
ultra-metal-poor halo stars [53] constrain the primordial 
metal abundance to be at least a few orders of magni- 
tude smaller than solar, we conclude that Planck will 
be unable to improve upon current constraints to the 
primordial metal abundance or, alternatively, that the 
standard CMB predictions for Planck are robust to pri- 
mordial metals at the levels allowed by current empirical 
constraints. 



II. EFFECT OF NEUTRAL METALS ON THE 
LYMAN-a DECAY RATE 



All metals (in the proper chemical sense of the term, 
i.e. not including noble gases, halogens and other non- 
metals) have a first ionization energy below 10.2 eV, 
which corresponds to the Lya transition in hydrogen. 
This means that neutral metals can provide continuum 
opacity in the vicinity of the Lya line by absorbing Lya 
photons in photoionization events. Since the photoe- 
jected electrons rapidly thermalize their energy, this re- 
sults in a net loss of resonant Lya photons, which would 
have otherwise been reabsorbed by ground state hydro- 
gen atoms. The presence of metals can therefore speed 
up hydrogen recombination by increasing the net rate of 
Lyman-a decays. A similar process was investigated for 
primordial helium recombination [24, 25, : in that case the 
presence of neutral hydrogen leads to continuum opacity 
in the He I 2^P° — 1^5* line. To estimate the impact 
of continuum opacity on the Lyman-a line, we use the 
analytic treatment presented for He I recombination in 
Ref. 1261. 



A. Continuum opacity in Lya due to 
photoionization of neutral metals 

The radiative-transfer equation in the vicinity of Lya 
for the photon occupation number including only true 
absorptions and emissions (i.e. neglecting resonant scat- 
terings) and continuum opacity is 



1 



Hvi^ya dt 



dv 



X2p 



= Tabs0(j^) , „ 



(1) 



where we approximated v « i/Lya in the prefactor on the 
left-hand-side. Tabs is the Sobolev optical depth for true 
absorption in the Lya line, (/>(i') is the line profile, and 
r]c is the continuum differential optical depth, given by 



_ nMOCCrpi(l^Lya) 



(2) 



In Eq. ([2]), n-j^o is the abundance of neutral metal M°, 
and (Tpi (p) is the photoionization cross section of M*^ at 
frequency u. We have assumed that (Jpi varies slowly over 
the Lya resonance (specifically, over the region which is 
optically thick for true absorption, which corresponds to 
a few tens of Doppler widths [H]), so we can approximate 
CTpi(i^) ~ crpi(z^LyQ)- Note that Eq. (fTl) assumes that the 
ionization state of M° is given by the Saha equilibrium 
equation (this translates in a ratio of continuum emission 
to absorption rates equal to / fv), even though 



this is not strictly correct (see Sec. II B) 

The net rate of 2p 
follows: 



Is decays is then obtained as 



( X2p 



\3a;i 



U]diy, (3) 



where the prefactor converts photon occupation numbers 
to photons per unit frequency per hydrogen atom, and we 
have approximated ly i^Lya in the multiplicative factor. 
Ref. [26j showed that the net decay rate in Lyman-a can 
be written in the following form: 



where 



std 



^""Ht^lya I X2p_ _ ^_„.,^„/T.. 



C-^7lH \2iX 



(4) 



(5) 



is the standard net decay rate in Lya, in the Sobolev 
approximation, for a large optical depth and assuming 
an incoming blackbody radiation field, and f is a correc- 
tion factor accounting for continuum absorption in the 
line. The correction factor £{tc) depends on the single 
parameter 



Tabsr2p?7c 

4772 



(6) 
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where T2p is the total inverse hfetime of the 2p state. 
The dimensionless parameter Tc can be interpreted as the 
continuum optical depth within the part of the Lya line 
which is optically thick for true absorption. For Tc — ?> 0, 
^{tc) 1, and for Tc > 0, £{tc) > 1, which is what one 
would expect as continuum opacity increases the net rate 
of decays in the line, as explained above. For <C 1, we 
have the following approximate expansion (see Eq. (117) 
of Ref. [Ml): 



f « 1 + 13tc 



< 1. 



(7) 



In order for primordial metals to change the net decay 
rate in Lya by ~ 1% (which is roughly the level de- 
tectable by Planck), we therefore need Tc ~ 0.001 near 
the peak of the visibility function. Extracting the rele- 
vant parameters from the multilevel atom code of Hirata 
|14j . we obtain, for z = 1100, 



« 0.7 X 10 



-3 0"pi('^Lya) 

10-i7cm2 



Xmo_ 
10-9 



(8) 



where Xm" = 't-m" /'^h is the abundance of neutral met- 
als relative to hydrogen. We see that for a characteristic 
photoionization cross section Cpi = 10"^^ cm^, a frac- 
tional abundance of neutral metals per hydrogen atom 
as low as ^ 10^^ would be potentially detectable. 



B. Ionization state of metals and results 

We now turn to the evaluation of the fraction of neutral 
metals /mo = timo/it-m- As a first approximation we use 
the Saha equilibrium value: 



(1-/mo) 



/mo 



5, 



Saha 



M 



ffM+ffe (27rmerr) 



3/2 



-Xm/T, 



. (9) 

where xm is the ionization energy of M", Uf, is the free 
electron abundance, and the g's are the degeneracy fac- 
tors for each species. For a standard recombination his- 
tory, at z = 1100, Eq. ([£]) gives (taking the ratio of degen- 
eracy factors to be unity) f^o = 5 x lO^'^, lO^**, 2 x 10"^ 
and 5 x 10"^ for xm = 10, 9, 8 and 7 eV, respectively, 
and we can anticipate that only metals with xm ^ 8 eV 
may have some impact on Lya. 

Saha equilibrium assumes that the ionizing radiation 
field is thermal. During hydrogen recombination, the ra- 
diation field develops large distortions in the vicinity of 
the Lya line, due to the slow escape of Lya photons (in 
fact, thermalization of these distorsions is so inefficient 
that they survive until today [37j). These non-thermal 
photons increase the ionization rate with respect to its 
thermal value, and the neutral fraction of metals is there- 
fore smaller than predicted by the Saha equation (see for 
example Ref. [9] for the case of lithium). The ioniza- 
tion state of the metal M is therefore rather determined 
by the balance of recombinations and photoionizations 
(this assumes the steady-state limit, valid so long as the 



photoionization rate is much larger than the Hubble ex- 
pansion rate, which is a very good approximation around 
the peak of the visibility function) 



(10) 



where aM is the — J-M*^ recombination coefficient and 
Pm is the M" ^-M^ photoionization rate. From Eq. (10 1, 
we obtain the neutral fraction, 



/mo 



/3: 



M 



(11) 



The photoionization rate I3m = Pcmb + /^dist comprises 
a thermal part, due to photoionizations by CMB pho- 
tons (from the ground state and excited states), which 
is related to the recombination coefficient through the 
detailed balance relation 



JCMB 



SuneOtM, 



(12) 



and a non-thermal part, due to photoionizations from the 
ground state by distortion photons. 



p(CMB) 



(13) 



where i^m = Xu/h and A/^ — fv — f^'^"^"' is the non- 
thermal distortion to the photon occupation number. 
The Lyman-a distortion peaks around z ^ 1400. We can 
therefore expect that distortions may start significantly 
affecting the ionization state of the metal M around red- 
shift z ~ 1400 xm/(10.2 eV). As a consequence, we ex- 
pect the Saha equilibrium approximation to be quite ac- 
curate around z ^ 1100 for metals with ionization thresh- 
old lower than ^ 8 eV. For the more interesting metals 
with Xm ^ 8 eV, however, spectral distortions will lower 
the neutral fraction with respect to the Saha value at 
z ^ 1100, making their detection more difficult through 
the effect considered here (we will consider the effect of 
additional free electrons due to the presence of ionized 



metals in Sec. IV I 



We have computed the neutral fraction of several met- 
als with atomic number Z < 26, using the fits of Ref. |28] 
for the photoionization cross sections, and the CHIANTI 
database for the recombination coefficients ^291 [3D] . We 
have extracted the Lya distortiorj^ from the two-photon 
code of Hirata [14] . We show the ionization state of beryl- 
lium (xbc = 9.32 eV), boron (xb — 8.30 eV) and silicon 
(xsi — 8.15 eV) as a function of redshift, for a standard 
recombination history, in Fig. [ij 



^ In principle, to be self-consistent, one should account for the 
continuum optical depth due to metal photoionization above Lyo 
and between i/m and vi^ya- Given that we find that this effect 
should not be detectable anyway, we have not implemented a 
more subtle treatment. 
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FIG. 1. Neutral fraction of beryllium, boron and silicon as a 
function of redshift. Thin lines represent the Saha equilibrium 
value given by Eq. ([9|. Thick lines represent a more accurate 
estimate accounting for distortions to the ambient blackbody 
field near Lyman- a. 
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FIG. 2. Minimum abundance of metals relative to hydrogen 
needed to result in a continuum optical depth in Lya Tc > 
0.001 at z = 1100. 



We show in Fig. [2] the minimal abundance of met- 
als detectable through its effect on Lya (i.e. such that 
Tc > 0.001 at redshift 1100). We see that the smallest 
detectable abundance would be xbc 3 x 10~^. Due 
to lack of data, we have not treated the case of other 
metals with xm > 9 eV, such as zinc (xzn — 9.39), ar- 
senic (xAs — 9.79) and gold (xau = 9.23), but do not 
expect significantly lower detectability thresholds unless 
they have unusually high photoionization cross sections. 



III. THE BOWEN 
RESONANCE-FLUORESCENCE MECHANISM 
FOR OXYGEN 



line. Neutral oxygen in the 2p^3d "^D^ state can then ei- 
ther directly decay back to the ground state or first decay 
to the 2p^3p state by emitting an infrared photon at 
1.13 /im, and subsequently cascade down to the ground 
state (in principle, atoms in the 2p'^3d ^D'^ state can also 
be excited or photoionized, but there are very few ther- 
mal photons energetic enough to do so) . The probability 
of the latter series of events (neglecting stimulated de- 
cays) is Pl.is^m = ^1.13Atm/(^1.13Ami + ^1025.76) ~ 0.3. 

Direct decays of excited oxygen to the ground state do 
not affect radiative transfer in the Ly/3 line, as they do 
not change the number of Ly/3 photons. On the other 
hand, absorptions in the 1025.76 A line followed by emis- 
sion of infrared photons degrade Ly/3 photons that would 
otherwise have been reabsorbed by neutral hydrogen. 
This effect is similar to the continuum opacity in Lya 
discussed in Section [llj except that this is now a reso- 
nant process. 

The escap^ probability in the Ly/3 line is enhanced 
by the probability that a Ly/3 photon is absorbed by O I 
(and then degraded) rather than by H L 



InoiA 



1025.76Pl.13Mm 



0.2 



3p,ls 



(14) 



where the multiplicative factors are the ratios of the de- 
generacy factors of the excited levels to those of the 
ground states. Note that the photon occupation num- 
ber redward of Ly/3 is slightly decreased by this pro- 
cess: /_ = a;3p/(3a;is) (1 — APcsc)- However, as long as 
APesc ^ 1 this has no detectable impact on radiative 
transfer redward of Ly/3. 

As oxygen and hydrogen have very similar ionization 
energies, we can assume that they have the same recom- 
bination history. More specifically, the ratio of ionized 
to neutral oxygen rapidly equilibrates to the correspond- 
ing ratio for hydrogen because the continuum above the 
ionization threshold is optically thick. We therefore have 
"-Oi/^-oii ~ |^iHi/"^Hii (the 9/8 comes from properly ac- 
couting for degeneracy factors, see Eq. (|9f) and as a re- 
sult we obtain noi/n-Hi ~ a;o/(l — Xe/'d) ~ xq- As long 
as xqi < 10~^ the damping wings of the O I 1025.76 
A line are optically thin, and the Bowen mechanism can 
only affect the net decay rate in the Doppler core of the 
Ly/3 line. The latter is very small anyway as the radia- 
tion field is very close to equilibrium with the 3p-ls ratio 
over many Doppler widths near line center. We have 
modified the escape probability from the Doppler core of 
Ly/3 in the recombination code HyRec i26J according to 
Eq. (14 1, and found that a minimal abundance of oxy- 
gen xq « 10^^ is required to affect the recombination 



There is an accidental coincidence between the hydro- 
gen Ly/3 fine at 1025.72 A and the O I 2p'^ ^P2-2p^3d ^D^ 
line at 1025.76 A. Ly/3 photons may therefore excite neu- 
tral oxygen instead of being reabsorbed in the hydrogen 



^ The term "escape" is somewhat misleading in this situation: pho- 
tons near the Ly/3 frequency do not actually escape more from 
the resonant region (in fact, their overall escape rate is even 
lower due to a slightly higher optical depth). They rather only 
"escape" reabsorption by neutral hydrogen. 
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history at a potentially detectable level A^e/a^e ~ 0.2% 
at z ~ 1100. Note that this would correspond to an en- 
hancement by a factor ^ 200 from the standard escape 
probability in the Ly/3 Doppler core. This stems from 
the fact that only a tiny fraction of Ly/3 decays ocurr in 
the Doppler core of the line, whereas most of them take 
place in the damping wings. The recombination history 
is therefore highly insensitive to the exact decay rate in 
the core. 



IV. ADDITIONAL FREE ELECTRONS DUE TO 
IONIZED METALS 

If metals remain ionized, they can contribute an addi- 
tional residual free-electron fraction at late times, Axg ~ 
xm+- In fact we have Ax^ = ■^xm+, as we show below. 
At late times the evolution of the free-electron fraction 
is given by 



-nuaBXeXp = -nuaBXeiXe 



xm+I 



(15) 



Eq. ( 15 ) is valid because the free electron fraction is many 



orders of magnitude above the Saha equilibrium value at 
late time (for a discussion, see Ref. [H]). If x^ is the 
unperturbed free electron fraction (i.e. obtained with 

xm+ = 0) and x^ — x^ + Axg, we obtain 



Aa;„ = 



a;° (2Aa;e - Xm+) , 



(16) 



which asymptotes to Axe — ^Xm+- The Planck satellite 
will be sensitive to fractional changes Axe/xe ^ 1% at 
late times. Since « 0.3 - 1 x 10"^ foj. 200 < z < 700, 
we conclude that a potential detection by Planck requires 
a fractional abundance of metals xm ^ 10~^ (in the case 
that metals remain fully ionized). Note that for a given 
fif,, the presence of metals also modifies the total abun- 
dance of hydrogen, nu, throughout the recombination 
history. However these modifications are degenerate with 
a mere change of fif, of Ine at the level of a few times 
10^^ and are therefore undetectable. 



V. CONCLUSIONS 

We have investigated whether a primordial metal con- 
tent could sufficiently affect the recombination history 
to be detectable in upcoming CMB data from Planck. 



We first considered the effect of photoionization of neu- 
tral metals by Lya photons. We showed that although a 
very small abundance of neutral metals would be enough 
to significantly affect the net decay rate in Lya, metals 
with ionization threshold below Lya are mostly ionized 
at z ^ 1100, and therefore undetectable. We also con- 
sidered the Bowen resonance-fiuorescence mechanism if 
primordial oxygen is present. This effect leads to an 
enhanced escape rate of Ly/3 photons and a speed up 
of recombination. We showed that it could lead to de- 
tectable changes for a primordial oxygen abundance of 
a couple hundredths of solar xq ^ 10^^. Finally, we 
pointed out that metals that stay ionized until late times 
provide additional free electrons and therefore change 
the late-time Thomson scattering optical depth. A frac- 
tional abundance a;M ^ 10~^ of primordial metals could 
be detectable through this effect. As a reference, the 
most abundant metal in the solar photosphere is oxygen 
{xo = 4.9 X lO-''), followed by carbon (xc = 2.7 x lO"*), 
neon (a;No = 8.5x10^^), nitrogen (xn = 6.8x10^^), mag- 
nesium {xug = 3.4 X 10~^), silicon {xsi — 3.2 x 10^^), 
iron (xpc = 3.2 x 10"^) and sulfur {xs = 1.3 x 10"''^). 
Other metals have fractional abundances xm < 10~^ in 
the Sun ^31) . As carbon, nitrogen, oxygen and neon are 
neutral at late times (due to their high ionization po- 
tential), we conclude that Planck could potentially de- 
tect primordial metals with an abundance at least a few 
tenths of solar. This is moreover an optimistic estimate, 
as the effect of metals is likely to be degenerate with the 
redshift of reionization or other cosmological parameters. 

Given that Lyman-alpha-forest measurements and 
ultra-metal-poor halo stars suggest a primordial metallic- 
ity much smaller than one hundredth solar, we conclude 
that the CMB can unfortunately not usefully constrain 
the abundance of primordial metals. At the same time, 
we also conclude that the CMB predictions for the Planck 
satellite are robust to a primordial metallicity allowed by 
current empirical constraints. 
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